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Abstract 

The azimuthal elliptic asymmetry v2 observed in heavy ion collisions, is usually 
associated with properties of the medium created in the final state. We compute the 
azimuthal asymmetry which is due to multiple interactions of partons at the initial 
stage of nuclear collisions, and which is also present in pA collisions. In our approach 
the main source of azimuthal asymmetry is the combination of parton multiple 
interactions with the steep variation of the nuclear density at the edge of nuclei. 
We apply the light-cone dipole formalism to compute the azimuthal asymmetry 
of prompt photons yield from parton- nucleus, proton- nucleus and nucleus- nucleus 
collisions at the RHIC energy. 
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1 Introduction 



Prompt photons, i.e. photons not from hadronic decays, are interesting since 
they do not participate in the strong interactions and therefore carry infor- 
mation about the initial state hard collisions. Nevertheless, measuring prompt 
photons is a challenge for experimentalists, partly due to the existence of 
large backgrounds coming from hadronic decays, which should be extracted. 
Even after this subtraction, there are other several sources for direct pho- 
tons, including thermal radiation from the hot medium and photons induced 
by final state interactions with the medium. In this paper we concentrate on 
the azimuthal asymmetry of prompt photons produced at the initial stage of 
relativistic nuclear collisions. 
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The PHENIX collaboration at RHIC has recently reported in Refs. [Tf2] the 
measurement of an azimuthal asymmetry of direct photon production, which 
has been studied recently in several theoretical papers [Bflfofo] . A novel mech- 
anism which produces an azimuthal asymmetry coming from the reaction's 
initial conditions was introduced in Ref . [3] . This is in contrast with the usual 
assumptions taken in approaches where the azimuthal asymmetry is only as- 
sociated with the properties of the medium created in the final state. We show 
that at least part of the direct photon azimuthal asymmetry, albeit small for 
AA and pA collisions, originates from initial hard scatterings between partons 
of the nuclei. In our approach, the main source of the azimuthal asymmetry 
originates from the sensitivity of parton multiple interactions to the steep vari- 
ation of the nuclear density at the edge of the nuclei, which correlates with 
the color dipole orientation. 

This paper is organized as follows: In sections 2 and 3 we introduce the main 
formalism for prompt photon production and discuss the relevance of color 
dipole orientation. In section 4, we introduce the azimuthal asymmetry for 
various collisions. In section 5, we present the numerical results. Some con- 
cluding remarks are given in section 6. 



2 Photon radiation in the color dipole formalism 



The transverse momentum (p?) distribution of photon bremsstrahlung, com- 
ing from the interaction of a quark with nuclear matter of thickness T^(&) = 
^^dzpA^z^b) (where the nuclear density pa is integrated along the parton 
trajectory at impact parameter b), integrated over the final quark transverse 
momentum can be written as [3~]T7] . 



d{lna)d l p T d l b {^) in jJ 

x F A (6,af!,ar 2 ,x), (1) 

where a denotes the fraction of the quark light-cone momentum carried by 
the photon, and 7g (a,r) is the light-cone amplitude for the qj fluctuation 
with transverse separation r. In this equation the QCD part is encoded in the 
function FA(b, afi, ar*2, x), which is a linear combination of qq dipole partial 
amplitudes on a nucleus at impact parameter b, 

F A (b, ar u af 2 , x) = Im/^(6, ar x ,x) + Im/^(6, af 2 , x) 

-Imf^ain -f 2 ),x), (2) 
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where the partial elastic amplitude can be written, in the eikonal form, in 
terms of the dipole elastic amplitude of a qq dipole colliding with a proton 

at impact parameter b [8], 



-jj d 2 s dm f£(s,r)T A (b + s)} 
l-exp[- f d 2 slmf q N Mf)T A (b + s}} 



(3) 

Dependence on the light-cone momentum fraction x of the target gluons and 
a are implicit in the above expression. The dipole partial elastic amplitude 
fqq was proposed in Ref. [3] to have the form 



= ^-/ ^^"^aSix, q, q') (e**> - e^ 1 ^) 
x (e^ - e-^^ 1 "^) , (4) 

where we defined a s = \Jct s (q 2 )a s (q' 2 ). The fractional light-cone momenta of 
the quark and antiquark are denoted by 77 and 1 — 77, respectively. The radiated 
photon takes away fraction a of the quark momentum. Therefore, we have the 
parameter 77 = 1/(2 — a) for the photon production. It is known that the center 
of gravity of qq is closer to fastest q or q. The generalized unintegrated gluon 
density_j T[x, q, q') is related to the diagonal one by T{x, q) = T[x, q,q = q'). 
Integrating over the vector s one can recover the dipole cross section cr^(r), 
and also r\ or a dependence will disappear 



2 j d 2 s!mf q N ^r) 
d 2 q 
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)a s (q 2 )J 7 (x,q). 



(5) 



Relying on the saturation shape of the dipole cross-section cr^(r) [TO] , the 
following form for a s T{x, g, g') was proposed in Ref. [3], 



a s F(x, q, q') = ^g^'2 jR 2 (x)e -i^(,)( 9 H 9 ' 2 ) e -|^(f-?') 2 ; (6) 

167T 

where the parameters <jq = 23.03 mb, Rq(x) = 0.4fm x (x/xo) 0144 with xq = 
3.04 x 10 -4 are fixed to HERA data for the proton structure function [10J. The 

1 This should not be mixed up with the generalized gluon density [9] which is 
off-diagonal in the longitudinal fractional momentum x. 
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parameter R 2 N = 5 GeV -2 is the t-slope of the pomeron-proton vertex [3j. The 
energy scale x which enters in the dipole amplitude is related to the measurable 
variable x = pr/w [llj . where w is the center of mass energy. Unfortunately, it 
is not possible to uniquely determine the unintegrated gluon density function 
from the available data. Nevertheless, the proposed form Eq. seems to 
be a natural generalization which preserves the saturation properties of the 
diagonal part [3j. After carrying out the integrations in Eq. (TjJ the dipole 
amplitude gets the following simple form, 



Im /il(s,f) 



'</</ 




[s + f(l -rj)]' 



2B e i 

[s + (1/2 - V )r\ 2 
2B P j 



+ exp 



[s — frjY 
2B P i 



(7) 



where we defined B e [ = R 2 N + Rq(x)/8. 



In Eq. ([T]), 7(? (a,r) is the light-cone (LC) distribution amplitude of the pro- 
jectile quark 7g fluctuation. Averaging over the initial quark polarizations and 
summing over all final polarization states of the quark and photon, we get 



inj 



07 9 («, n)(f) jq (a, r 2 ) = -^mJa 2 {a 2 i^o(amgri)K (am g r 2 ) 



+ [1 + (1 - aY^—K^am^K^am^)}, (8) 
r x r 2 



where K 0t i(x) denotes a modified Bessel function of the second kind and m q 
is an effective quark mass, which can be conceived as a cutoff regularization. 
We take m q = 0.2 GeV for the case of direct photon production [TT]. It has 
been also shown that a value of m q = 0.2 GeV is needed in order to describe 
nuclear shadowing effects [12J. 

Expression ([1]), with the exponentials expanded to first order in the nuclear 
thickness, provides also the cross-section for direct photon production in hadron- 
hadron collision. We have recently shown that in this framework one can obtain 
a good description of the cross section for prompt photon production data for 
proton-proton (pp) collisions at RHIC and Tevatron energies [llj. Predictions 
for the LHC in the same framework are given in Ref. [15] ; while to compare 
with the predictions of other approaches at the LHC see Ref. [T4] . 
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Fig. 1. The partial elastic amplitude Ini/A (mb) of the qq dipole on a proton at 
impact parameter s as a function of dipole size r and angle (5 defined between s and 
r for two values of a = 0, 1. We use a fixed value of x = 0.01 for all plots. 

3 Azimuthal asymmetry and dipole orientation 

The main source of azimuthal asymmetry in the amplitude is the inter- 
play between multiple rescatering and the shape of the physical system. The 
key function which describes the effect of multiple interactions is the eikonal 
exponential in Eq. ([3]), while the information about the shape of the system is 
incorporated through a convolution of the impact parameter dependent partial 
elastic amplitude and the nuclear thickness function. Notice that the initial 
space-time asymmetry gets translated into a momentum space anisotropy by 
the double Fourier transform in Eq. ([1]). 

It is quite obvious intuitively that although the Rutherford scattering cross 
section is azimuthally symmetric, the azimuthal angle of the radiated photons 
transverse momentum at a given impact parameter s correlates with the di- 
rection of s. In terms of the partial elastic amplitude f^(s,r), it means that 
the vectors f and s are correlated. This is the key observation which leads to 
an azimuthal asymmetry in pA and AA collisions. 

In Fig. we show the partial dipole amplitude f^(s,r) as a function of the 
dipole size r and the angle (3 between s and r, at various fixed values of s for 
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two values of a = 0,1. Notice that the generic feature of the partial dipole 
amplitude, e. g. its maximum and minimum pattern, changes with a. One can 
see that for very small dipole size r the dipole orientation is not important. 
For very large dipole size r compared to the impact parameter s or very small 
values of s the dipole orientation is also not present. Notice, however, that it is 
not obvious a priori how the convolution between the partial dipole amplitude 
and the nuclear profile, which leads to a more complicated angle mixing, gives 
rise to a final azimuthal asymmetry. The main aim of this paper is to calculate 
such azimuthal asymmetry without using any approximations. 



4 Azimuthal asymmetry in qA, pA and AA collisions 



The azimuthal asymmetry of prompt photon production, resulting from parton- 
nucleus (qA) collisions, is defined as the second order Fourier coefficients in 
a Fourier expansion of the azimuthal dependence of a single-particle spectra 
Eq. ([!]) around the beam direction, 



where the angle is defined with respect to the reaction plane. Some of inte- 
grals in the above expression can be analytically performed. After some tedious 
but straightforward calculation one obtains 



where the functions typ and which contain information about ^q fluctu- 
ation, are defined by 



/^#cos(20) 



da" A (q^qj) 




(9) 




f£° drrfy D (p T ,r,a)$ D (b, r,a) + 2irAf(a,p T ) 



(10) 
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, T , / n a em 2 4 / J 2 (p t r)K (am q r) r 

+ [1 + (1 - a) 2 ] f 2 , 2 (^iM - UPtr)) Ki(am q r) 

TOlTfl \ 

+ J 2 {ptr)K (am q r) ^-J 2 (p t r)Ki(am q r) J }; 

, T , / \ «em ro 2 4 / J {p t r) K {am q r) r \ 

+ J {ptr)K {am q r) ^ J (p t r)i^i(am g r) J }; 



27r 2 ^(p^ + (am,) 2 ) 2 (pf, + (am,) 



2^2 



(11) 



where J n (x) , n = — 3 denotes the Bessel functions of the first kind. The 
functions $jv,i? (b, r, a) in Eq. f TTOT) contain information about the QCD part 
and also the shape of the nucleus; 



$ N (b, r,a) = - d(3e-I d2g W£(s>-)t a (M cos ( 2 ^) ; 

J — 7T 

$ D (6, r, a) = r df3e~I d2g W^>^t a (M ; ( 12 ) 

where /3 is the angle between b and r, and the dipole amplitude /^(s, r, ) was 
defined in Eq. ([TJ). It is interesting to note that in the final form of v 2 A , Eq. ffTOj) . 
the angle between the impact parameter b and the transverse momentum 
of the projectile quark p? disappeared, and then the azimuthal asymmetry is 
directly related to the dipole orientation with respect to the impact parameter 
b through the angle (3 (see Eq. ( TT2]) ). Therefore, if one neglects the dipole 
orientation the azimuthal asymmetry becomes identically zero, regardless of 
both the nuclear profile and the dipole cross-section parametrization. 

In order to obtain the hadronic cross section from the elementary partonic 
cross section Eq. ([TJ), we use the standard convolution based on QCD factor- 
ization [T5] . 



da^(pA -> iX) 1 f 1 da x 1 da qA {q -> qj) . 

- F 2{ — ,Pt)— — : — -=i (13) 



dx F d?p T d?b x 1 +x 2 J Xl a a d(lna)d?p T d?b 

We take the parametrization for the proton structure function F 2 (x, Q 2 ) given 
in Ref. [TEj. Here x\ denotes the fraction of the light-cone momentum of the 
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projectile hadron carried away by the photon, and we define %2 = x\ — xf, 
where xp = 2pi/^/s is the Feynman variable. The azimuthal asymmetry of 
photon yield in proton-nucleus (pA) interactions is then defined as 



v 2 



(b,Pr) 



dtr~i(pA^-yX) 
dxpd?pTd?b 



(14) 



The above equation can be also simplified to an expression similar to Eq. (I10p . 
but now augmented with the proton structure function and an extra integral 
over the variable a. 

The spectra of photon bremsstrahlung from nucleus-nucleus (AA) collisions 
can be obtained from Eq. ffl3|) by weighting the cross-section with the density 
overlap factor of nuclei. In principle, one may obtain the cross-section for AA 
collisions in the same fashion as was done in Eq. ( fl3l) for the case of pA 
collisions, that is by making a convolution with the nucleus structure function 
instead of the proton structure function. However, the medium modification 
of nucleon structure function for the range of pt values we are interested in 
is less than 20% [18] and therefore will not change the overall prediction. The 
azimuthal asymmetry of photon yield from collisions of two nucleus Ai and 
A 2 at impact parameter B is defined as 



vt lA2 (B,p T )- 
Gn : 



/^0cos(20)£ 



N 



d<p G D 1 

d o cos(29i) =; — T A2 (O2) 

dxpd 2 pTd 2 b\ 

+ / crocos(29 2 ) — T Al {bx); 



dxpd 2 pTd 2 b2 



d 2 b 



'da^(pA l -> 1 X) 
, dxFd 2 prd 2 bi 



T A2 (h) + 



da\pA 2 -f n,x) 
dxFd 2 pTd 2 b2 



(15) 



where we used the notation 62 = b + B, bi = b (b is the impact parameter of 
the pAi collision) and the angle 0i (62) is the angle between the vectors 61 ( 
b 2 ) and B, respectively. The factor cos(29i j2 ) in the above equation relates 
the reaction planes of the pA and the AA collisions. The integral over b in 
Eq. (!T5]) covers the almond shape area of the nucleus- nucleus overlap. 
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Fig. 2. The azimuthal anisotropy of prompt photon production coming from 
quark-nucleus collisions, at various impact parameter b and at the RHIC energy, for 
both the Woods-Saxon (WS) and hard sphere (HS) profiles. The relative fraction 
of the quark momentum carried by photon is taken to be a = 1 for all curves. 



5 Numerical result and discussions 



We will perform a numerical calculation for the RHIC energy y/s = 200 GeV, 
at midrapidities. The only external input is the nuclear profile. First, we take 
a popular Woods-Saxon (WS) profile, with a nuclear radius Ra = 6.5 fm and 
a surface thickness £ = 0.54 fm, for Pb+Pb collisions |17j . 

In Fig. (j2J) we show the calculated values of v\ defined in Eq. (P), for fixed 
a = 1, at various qA collision impact parameters b, and at the RHIC energy. 
For central collisions, the correlation between nuclear profile and dipole orien- 
tation is minimal. In fact, if the nuclear profile function was constant, then the 
convolution between the nuclear profile and the dipole orientation, defined in 
Eq. (J3j), would be trivial, and v% A becomes then identically zero. Therefore, the 
main source of azimuthal anisotropy is not present for central collisions. This 
can be seen in Fig. (12"]) , where a pronounced elliptic anisotropy is observed for 
collisions with impact parameters close to the nuclear radius Ra, where the 
nuclear profile undergoes rapid changes. Therefore, an important parameter 
which controls the elliptic asymmetry in this mechanism is \b — Ra\- We have 
verified this numerically by taking different Ra values for the WS profile, but 
with the same surface thickness £. 

In Fig. ([3]) we show v% A for prompt photons produced in non-central pA col- 
lisions, for the WS nuclear profile. Notice that at very small transverse mo- 
mentum our results are not reliable, because the employed proton structure 
function is not valid. 
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Fig. 3. The impact parameter dependence of prompt photon azimuthal asymmetry 
for proton-nucleus collisions at various impact parameter b at RHIC energy. For the 
nuclear profile, we have taken the Woods-Saxon (WS) profile for all curves. 

In Fig. we show the prompt photon azimuthal asymmetry for AA collisions 
defined in Eq. (fT5l) . at various impact parameters and at RHIC energy. The 
absolute value of v AA turns out to be reduced compared to both v% A and 
The reason is that the integrand in Eq. (115p gets contributions only from 
semi-peripheral pA collisions where our mechanism is at work, and most of 
the integral over b does not contribute. This significantly dilutes the signal. 
For prompt photon, there is no suppression mechanism related to medium 
effects, as for the case of hadron production in central AA collisions compared 
to pp collisions [32]. For hadronic t> 2 , the inclusion of such suppression effects 
significantly enhances the azimuthal asymmetry coming from this mechanism 
[2U] . The other diluting factor for v^ A is the presence of an extra cos(2G) in 
Eq. (fl~5|) . which accounts for the changing of the reaction plane going from pA 
to AA collisions. At high p?, where the dipole size is very small, the dipole 
orientation becomes less important and consequently the correlation between 
the dipole cross-section and the nuclear profile disappears, i.e. the azimuthal 
asymmetry vanishes. This can be seen from Figs. (|2"]I5P| . where t>2 for all qA, 
pA and A A collisions approaches to zero at high p?. 

In our approach the profile of nuclear density at the edge is a very important 
input, since the elliptic asymmetry stems from the rapid change of nuclear 
density at the edge. In order to show this more clearly and to estimate the 
theoretical uncertainty of our calculations, we obtain the elliptic asymmetry 
for the hard sphere profile with a constant density distribution, pA = PqQ{Ra~ 
r), with the same nuclear radius Ra as the WS profile for Pb. In Fig. ([2]) we 
show the azimuthal asymmetry v% A for the hard sphere profile (HS) at impact 
parameter 6 = 7 fm. For the HS profile the nuclear thickness changes at nuclear 
radius more steeply compared to the WS profile and consequently the elliptic 
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Fig. 4. The impact parameter B dependence of prompt photon azimuthal asymme- 
try, for Pb+Pb collisions at RHIC energy, with the WS nuclear profile. 

asymmetry is significantly bigger. Notice that at impact parameter b = 7 
fm for the HS profile, even though there is no matter, the correlation between 
color dipole and the variation of nuclear thickness still exists and consequently 
the azimuthal asymmetry is not zero. 



6 Summary and final remarks 

We have computed the azimuthal asymmetry of prompt photons originating 
from primary hard scatterings between partons. This can be accounted for by 
the inclusion of the color dipole orientation, which is sensitive to the rapid 
variation of the nuclear profile. We showed that the azimuthal asymmetry t> 2 
coming from this mechanism changes the sign and becomes negative for pe- 
ripheral collisions, albeit it is extremely small. The first experimental attempts 
by PHENIX to extract the elliptic flow of direct photons yielded results which 
are compatible with zero within error bars pp. However, the data was still 
contaminated with background from hadron decays. Recently, the PHENIX 
collaboration has presented preliminary data on v 2 of direct photons [2J. Al- 
though the systematics errors are still very large, the data indicates that v 2 of 
direct photons is larger than our prediction for v 2 of prompt photons. There- 
fore, if data is confirmed, other sources for v 2 of direct photons should be also 
important. 

There are, however, a number of caveats in our approach which need further 
study before taking the numbers predicted here at face value. As we have 
shown, the tail of nuclear profile is an important external input in this mech- 
anism and quite significantly affects the results. For example, we showed that 
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the azimuthal asymmetry is enhanced almost by an order of magnitude for 
the HS nuclear profile compared to the WS one (see Fig. <^). Unfortunately 
the tail of all available nuclear profile parametrizations is less reliable, since 
it is not well probed by electron scattering and is obtained by a simple ex- 
trapolation. This is also due to the fact that the neutron distribution, which 
may be more important on the periphery, cannot be properly accounted for by 
electron scattering data. This brings uncertainty to our results. Furthermore, 
as we already pointed out, due to lack of experimental data, there is some 
freedom left to define the off-diagonal part of the unintegrated gluon density. 
Although it is not important for the total cross-section, it plays an important 
role for the azimuthal asymmetry. 
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